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Abstract 

Both direct and indirect weak nonresonant interactions are shown to produce entanglement 
between two initially disentangled systems prepared as a tensor product of thermal states, provided 
the initial temperature is sufficiently low. Entanglement is determined by the Peres-Horodeckii 
criterion, which establishes that a composite state is entangled if its partial transpose is not positive. 
If the initial temperature of the thermal states is higher than an upper critical value T uc the minimal 
eigenvalue of the partially transposed density matrix of the composite state remains positive in 
the course of the evolution. If the initial temperature of the thermal states is lower than a lower 
critical value T\ c < T uc the minimal eigenvalue of the partially transposed density matrix of the 
composite state becomes negative which means that entanglement develops. We calculate the lower 
bound To, for T] c and show that the negativity of the composite state is negligibly small in the 
interval Tn, < T < T uc . Therefore the lower bound temperature To, can be considered as the critical 
temperature for the generation of entanglement. 

PACS numbers: 03.67.Mn,03.65.Ud,03.67.-a 
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I. INTRODUCTION 



Efficient simulation of quantum dynamics on classical computers is hampered by the 
problem of scaling: the complexity of computation in quantum dynamics scales exponentially 
with the number of degrees of freedom |l[. The reason for this exponential growth is the 
entanglement of the degrees of freedom that is generated during the evolution. This problem 
is of a fundamental character: entanglement is viewed as one of the main peculiarities of 
the quantum dynamics as compared to its classical counterpart . Asking under what 
conditions entanglement is generated along the evolution of the quantum system is closely 



associated with the question of the quantum-classical transition 



It is customary in quantum dynamical simulations to assume that the initial state of 
the composite system is factorized in the relevant local basis 6]. An important question is 
whether the product form is conserved along the evolution j?[ 8|. The answer was generally 



found to be negative both for the pure 



and for the mixed state dynamics. 
A pure composite state is entangled if and only if it is not factorized in the local basis. 
For mixed states the situation is more complex Q|. For a bipartite composite system 
separability 11] is defined as a decomposition of the density matrix of the composite system 
in the form 

P12 = ^PiPi ® P2> (!) 

where < Pi < 1 and YliPi = 1 anc ^ Pi an d P2 are density matrices on Hilbert spaces 
of the first and the second subsystem, respectively. Separable states exhibit only classical 
correlations. States that cannot be represented in the form exhibit correlations that 
cannot be explained within any classical theory and are said to be entangled. There are 
two qualitatively different kinds of the mixed states entanglement [l^ : free entanglement 
and bound entanglement. Free entanglement can be brought in a form useful for quantum 
information processing and bound entanglement is "useless" in this sense. 

Separable states are not of the product form generally. Thus the important question 
remains, under what conditions the mixed state of the composite system evolving from the 
initial product (or generally separable ) state develops entanglement along the evolution. If 
quantum correlations in the composite system do not develop during the evolution one may 
speculate that the dynamics of the composite system has classical character. A possible 
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practical implication is that this "separable dynamics" could be simulated efficiently on 
classical computers. 

The dynamics of entanglement was investigated recently in various systems: the quantum 
Brownian particle harmonic chain two-qubits system interacting with the common 



19, 



harmonic bath [ljj], Jaynes-Cummings model lid], NMR various spin systems 18, 
I20I ]. Morse oscillator coupled to the spin bath |2JJ and bipartite Gaussian states in quantum 
optics (2^] to mention just some cases. The purpose of the present paper is to explore the 
temperature dependence of entanglement generation in the course of evolution of a bipartite 
state in the limit of weak coupling and nonresonant interaction between the parts. Under 
these limitations nondegenerate perturbation theory was applied to the calculation of the 
bipartite entanglement in the evolving composite system. We have considered two cases 
of interaction - (1) direct interaction, when two initially disentangled systems are brought 
into contact ( Cf. FigGJ, and (2) indirect interaction, when two noninteracting and initially 
disentangled systems are brought into contact with the third party (Cf. FigEJ). In each case 
the initial state of the composite system was taken to be the product of the thermal states 
of the parts. 
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To establish quantum entanglement the Peres-Horodeckii criterion is employed 
The Peres-Horodeckii criterion states that the bipartite system is entangled when the par- 
tially transposed density matrix of the system possesses a negative eigenvalue. The converse 
statement is generally not true: there exist inseparable states whose partially transposes 
are positive It is proved in Ref.|3| that states whose partially transposes are positive 
(PPT states in what follows) do not exhibit free entanglement. Therefore PPT states are 
either separable or bound entangled and as a consequence are not useful in quantum infor- 
mation processing. In the context of simulating a quantum composite system with classical 
computers, we are interested in the possibility of maintaining a separable form (Cf. Eq.((TJ)) 
during the evolution. We conjecture on the basis of Ref. 26], where it is proved that PPT 
density matrices of sufficiently small rank are separable, that for a state that remains PPT 
during the evolution separability can be obtained by embedding in a larger Hilbert space. 

Applying the Peres-Horodeckii criterion to the case (1) we show that for sufficiently low 
initial temperature of the subsystems the interaction does induce entanglement unless the 
ground state of either one of the subsystem is invariant under the interaction. A lower 
and upper critical temperatures T} c and T uc exist such that if the composite system evolves 
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from the initial thermal state with temperature T < T\ c the minimal eigenvalue of the 
partially transposed density matrix becomes negative in the course of the evolution and if 
T > T uc the minimal eigenvalue of the partially transposed density matrix stays positive. 
The lower bound TJ& of the lower critical temperatures TJ C was calculated in the limit of 
weak intersystem coupling and shown to be tight: the negativity of the composite state 
(27 1 . which is a quantitative counterpart of the Peres-Horodeckii criterion and a measure of 
entanglement, is shown to be generally negligible for temperatures in the interval T«, < T < 
T uc . Therefore, according to the Peres-Horodeckii criterion, when T < 7}& the composite 
system develops entanglement in the course of the evolution and when T > the composite 
state remains PPT state. 

The question addressed in case (2) of indirect coupling is what are the conditions on the 
interaction with the common bath and on the initial temperature of the states which cause 
entanglement of the noninteracting systems? Two scenarios with time scales separation are 
studied: (a) two "slow" noninteracting systems coupled to a "fast" third party (b) two "fast" 
noninteracting systems coupled to "slow" third party. Under some technical assumptions 
about the form of the interaction we find in both cases that for sufficiently low initial 
temperature of the noninteracting systems entanglement is induced by the interaction with 
the third party. We calculate the lower bound temperature To, in both cases of the time 
scales separation. In the system of two noninteracting spins, coupled to the common bath, 
the lower bound coincides with the T uc . 

In both cases (1) and (2) the evolution starts from an uncorrelated initial state of the 
composite system represented by the tensor product of thermal states of the subsystems in- 
volved. As a consequence, initially the eigenstates of the partially transposed density matrix 
of the composite state are nonnegative. The evolution under the interaction perturbs the 
initial state. The new eigenvalues of the partially transposed density matrix are calculated 
by the nondegenerate perturbation theory assuming the coupling is weak and the interaction 
is nonresonanant. The time dependence of the minimal eigenvalue is not analyzed in detail. 
As the time evolution of the density matrix is quasiperiodic the minimal eigenvalue of the 
partially transposed density matrix is also a quasiperiodic function of time. The interaction 
is said to induce entanglement if the minimal eigenvalue becomes negative in the course of 
the evolution. 
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FIG. 1: The coupling scheme for two directly interacting systems. 



II. ENTANGLEMENT BETWEEN TWO DIRECTLY INTERACTING SYSTEMS 



A composite system A <g> B evolves under the following Hamiltonian : 



H 



total 



H + 7 V. 



(2) 



where H = H a 1 + 1® V = V a <g> Y b and 7 scales the magnitude of the interaction. 
Let the initial state be 



P(0) 



Pa ® Pbi 



(3) 



where both p a and p b are thermal states: p ab = Z~\ exp(— H^/T), where Z~\ is the nor- 
malization factor. The Boltzman constant ks is unity throughout the paper. The evolution 
is followed in the interaction picture. Then 



~dt 



Va(t)®V b (t),0 



(4) 



where 



p'(t) = e~ mt p e iH ', 



(5) 



Here and in the rest of the paper we take ft = 1. It is clear that the density matrix pit) 
is separable if and only if p{t) separable. In what follows the tags in p{t) are omitted for 
simplicity. In the first order in the coupling 7 the evolution of p becomes: 

p[t) = 0(0) - i! [ \v a (t') ® y b (t'), p(0)l dt'. (6) 

Entanglment of pit) is established by the application of the Peres-Horodeckii criterion. 
This is carried out by calculating the partial transpose of the state. The partial transposition 
T a with respect to subsystem A of a bipartite state p ab expanded in a local orthonormal 
basis as p ab = p%j,u K) (j\ ® \k) {l\ is defined as: 

The spectrum of the partially transposed density matrix does not depend on the choice of 
local basis or on the choice of the subsystem with respect to which the partial transposition 
is performed. By the Peres-Horodeckii criterion the eigenvalues of a partially transposed 
separable bipartite state are nonnegative. 

The density operator (jHJ) under the partial transposition (T a ) becomes: 

pit) Ta = p(0) Ta -*l I \v a (t')®V b (t'),p(0)} Ta dt'. (8) 



Let \ik) = \i) <8> \k) be the local orthonormal basis of the system A £g> B composed of the 
eigenvectors of the unperturbed Hamiltonian H = H a + H&: 

H a , fe |<> = El b |<) , (9) 

where E l a b , i = 1,2,..., is the unperturbed energy spectrum of the Hamiltonian H 0j fe. The 
initial state is of the tensor product form, Cf. Eq.(JHJ), therefore : 

p(0) T " \ik) = p(0) \ik) = P a ®P b \ik) = P lk \ik) , (10) 

where P ik = p a ,iP b ,k and p a ,i,Pb,k are defined by p a>i = (i\p a \i) and p b)k = (k\p b \k). The 
matrix elements of p{t) Ta in the chosen basis are given by: 



(ik \p{t) Ta j jl) = P ik S(ik),(ji) + Mucjt, (11) 

where 
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M, 



ik,jl 



/j I ( ik 

i'J J ((ik 



v a (f)®v b (t'),p a ®p b ) T ° ji)dH 



p a V a (t') T ®V b (t')p b 



jl 



V a (t') T p a ®p b V b (t') 3 l))dt' 



— 27 (ik 

= ij{Pu-Pjk) J U 
= 7 (P« - P ife ) (j 



(12) 



V,. 



v a (t') T 

'jfe 



J )(k 



l)dt' 



V, 



V, 



it(El-El+E' b -E*) _ 1 



A- 

/ (E* - £S + ^ - £*) 



where X designates the transpose of the operator X. 

When T = 0, the zero eigenvalue of the initial state p(0) is degenerate. As a result the 
zero eigenvalue of the partially transposed initial density operator p(0) Ta = p(0) is also 
degenerate. The zero eigenvalues correspond to empty initially unoccupied states. By the 
standard secular perturbation theory the first order correction to the degenerate eigenvalue 
A^ = of the matrix p(0) Ta is given by 



\M nn i — X^5 nn i I 



0. 



(13) 



where \n) and \n') are eigenvectors of the matrix pj ® p b} corresponding to the degenerate 
A(°) = 0. Since p(0) Ta = p^ ® p h — p a ® p 6 the eigenvectors of pj ® p 6 , corresponding to 
A< > = are |n) = {|1) ® |i) , |z) <g> |1) , |z) ® [?') ,i, j = 2,3, ..}. 

Therefore at T = 0, P^ = S^Ski, and by inspection of Eq. (|T2"j) . the only nonvanishing 
matrix elements in the degenerate subspace spanned by \n) and \n') are Mijji and Mji t u 
where either i ^ 1 or j ^ 1. Since the trace of the matrix M is zero, either all its eigenvalues 
vanish or some of them are negative. All the eigenvalues of M cannot vanish unless M — 0, 



which from Eq. (fT2"|) implies V a , p a 



or 



V b ,p 6 



0, i.e. the ground state of either 



one of the subsystem is invariant under the interaction. In this case the interaction acts 
locally on the subsystems and cannot entangle them. Otherwise there are negative solutions 
to Eq. (fT3*|) and as a consequence the partial transpose of the density operator attains 
negative eigenvalues already in the first order in the coupling. Therefore, according to the 
Peres-Horodeckii criterion, entanglement develops at zero temperature. 



To simplify the study of the generation of entanglement at finite temperatures it is as- 
sumed that the only non zero matrix elements of V aj & are those between neighboring states, 
i.e. (V a ,b)ij °^ 8ij±i- Under this assumption the partially transposed density matrix p(t) Ta 
obtains the following structure: 

/ r, n n m. , „„ n n n \ 



P(t) 



To 



Pu 








M llt22 














Pl2 


Ml2,21 





-^12,23 











^1*2,21 


P21 








-^21,32 





/If* 

iw ll,22 








-P22 








-^22,33 





/If* 

1V1 12,23 








-^23 


^23,32 











/If* 

iw 21,32 





/If* 

iw 23,32 


P32 














Af* 

JW 22,33 








-^33 



(14) 



v 

where Pij are defined after the Eq- ffTUj) and Mkiji by Eq. (fT2]) . 

There are two kinds of matrix elements Mkiji' ^fci,(fc+i)(i+i) and M hi ^ k ±i)(i^i) (other 
elements are their counterparts under the transposition). Matrix elements M k i t (k+i)(i+i) 
couple the unperturbed eigenvalues Pki and P(k+i)(i+i)- For small coupling strength 
7 |-^fci,(fc+i)(i+i) I ^ Pki and the contribution of M ki> ^+1)^+1) to the correction to Pki 
is negligible and cannot make the eigenvalue negative. On the other hand, the ratio 
\M ki ,{k+i){i+\)\ /P(k+i)(i+i) «: 7(P fc (i+i) - P{k+i)i) / P(k+i){i+i) can in general be arbitrary large 
for low temperatures but for sufficiently high temperatures it tends to zero and as a conse- 
quence the contribution of the coupling element Mki^k+i)(i+i) to the correction to P^+i)(i+i) 
is negligible. It can be checked along the same lines that the ratio of the coupling matrix 
elements M fci) ( fc±1 )(j T i) to the unperturbed eigenvalues Pki and P(k±i)(i^i) of the partially 
transposed density matrix (fT4^) vanish for sufficiently high temperature. Therefore, at least 
for composite systems with finite Hilbert space dimensions, there exists a finite upper crit- 
ical temperature T uc . Above T uc the spectrum of the partially transposed density matrix 
remains positive (PPT). In close vicinity of T uc from below the minimal eigenvalue becomes 
negative in the course of the evolution. These conclusions stay in accord with a general 
result [3, that finite dimensional composite states in sufficiently small neighbourhood 
of the maximally mixed state (i.e. thermal states at infinite temperature) are separable. We 
conjecture, that for an infinite composite system, the upper critical temperature T uc exists 
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if the energy spacing is bound. 

At sufficiently low initial temperature the minimal eigenvalue of the partially transposed 
density matrix becomes negative in the course of the evolution. This means that there exists 
a finite lower critical temperature T) c . Below T) c the composite systems A <g> B develops 
entanglement. In sufficiently close vicinity of TJ C from above the state remains PPT in 
the course of evolution. It is possible that T\ c = T uc . This equality is confirmed in all 
numerical tests. A lower bound 7]& for the lower critical temperature T\ c can be calculated 
using perturbation analysis. It is shown that this bound is tight since the free entanglement 
in the interval < T < T uc is negligibly small under the weak coupling assumption. 
Therefore, from the practical point of view the lower bound T/& for 2] c can be considered 
as the critical temperature for entanglement. For simplicity the lower bound for the lower 
critical temperature is termed "the lower bound temperature" throughout the paper. 

At low temperatures the leading order contribution to the negative eigenvalue of the 
partially transposed density matrix comes from the matrix elements Mu^, Mi2,2i (and 
their complex conjugates) that do not vanish at T = 0. Therefore, to the leading order in 
7, the nonvanishing eigenvalues of the partially transposed density matrix Eq. (jl4j) are the 
eigenvalues of the following effective partially transposed density matrix : 

/ Pn M llj22 \ 

Pl2 Mi2, 2 l 

M* 2>21 P 21 
\M* lj22 P 22 / 
The critical temperature, calculated for the effective 4x4 matrix (|15|). is a lower bound for 
the lower critical temperature T; c of the bipartite system A <g> B. The eigenvalues of Eq. 
()15|) are eigenvalues of two 2x2 matrices: 

P\2 Ml2,21 
M*2 2l Pl\ 



p{t) Ta 



eff 



(15) 



(16) 



and 

.22 

^11,22 ^22 

The eigenvalues of the matrix (jlfi)) are: 



(17) 



P 12 + P 21 ± y/(P 12 + P 21 ) 2 - 4(P 12 P 21 - |M 12 . 21 | 2 ) 

A ± = o ' y LH ) 



where from Eq.()12j): 

M X2M 



7(2 



V, 



V, 



p itAEn _ i 



P 



22) 



(19) 



where we define AEu = E\ — E\ + El — E^, which is the lowest joint excitation energy of 
the composite system. 

From Eq. (fTHj) . A_ will be negative whenever P12-P21 < |^i2,2i| 2 an d positive if 
-P12-P21 > |^i2,2i | 2 - The lower bound temperature T[ b is evaluated from the condition 
-P12-P21 = |Mi2, 21 1 2 - Since |Mi 2 ,2i| is an oscillating function of time (Cf. Eq. (fT9|) ) the 



amplitude of |M 12i 2i| is taken to be equal to y/PwPii- 
2 7 



V, 



(Pi 



AE U 

Assuming that T ih is low P n — P 22 ~ P\\ and then 

2 7 



11 



P 



22 1 



P\iPz\- 



(20) 



AE 



11 



V, 



V, 



' P 2 P 3 



(21) 



Since e isa monotonic function of the temperature, at T > A_ > and at T < 
A_ < 0. Finally, the expression for the lower bound temperature becomes: 

AE n 



Tu 



v. 



(22) 



So far only two of the eigenvalues of the matrix (|15|) have been evaluated. The other two 
eigenvalues are found to be strictly positive at and above the temperature Tu,. Therefore, 
the expression (}2*2*|) defines the critical temperature for the partially transposed effective 
density matrix (|T5|) and the lower bound temperature of the partially transposed density 
matrix ()14j) . 

Eq. (|2*2*j) can be generalized to an interaction term of the form ^ 7iV a ® V b : 



T, 



AE U 



ib 



2 In 



AEi 



E,7,(2 



V, 



(23) 



provided Ylili (2 V a -0 (2 V b l\ 7^ 0. When this term vanishes there is no entanglement 
in the first order in the coupling strength 7. 

For the system of two interacting spins the lower bound Tu, given by Eq. (}2*2*j) coincides 
with the upper critical temperature T uc therefore in this case the critical temperature exists 
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FIG. 2: The shaded area in the parameter space of the inverse initial temperature T of two spins and 
the logarithm of the inverse coupling strength 7, represents values of T and 7, where entanglement 
does not develop in the course of the evolution. The composite system of two spins evolves from 
the initial product of thermal states under the Hamiltonian H = ^u;(<x" <g> i + (\/2 — l)i <g> a b z ) + 
7(6-" (&a h x — &y(& a h y ). The evolution is calculated numerically for u = 1. The border of the shaded 
area represents T uc calculated numerically. The dashed line represents To, according to Eq. (|22|). 
Up to the coupling 7 = 0.1 approximates TJ C very well. 

in the strict sense. Fig. 121 shows results of the numerical calculation of the critical tem- 
perature as a function of coupling strength for a system of two interacting spins evolving 
from the initial product state of two thermal states. The Peres-Horodeckii criterion was 
used and the partial transpose of the evolving density matrix was calculated numerically 
to determine entanglement. The shaded area in the parametric space of the logarithm of 
inverse coupling and the inverse initial temperature represents the values of the parameters 
where no entanglement develops. For coupling up to 7 = 0.1 given by Eq. (J22j) (the 
dashed line) corresponds well to the numerical values of T uc . It is interesting to note that 
for large values of coupling the critical temperature asymptotically tends to a finite constant 
value of the same order of magnitude as the energy difference between the first excited and 
the ground state of the unperturbed composite system. 

At T < Tib the minimal eigenvalue of the partially transposed state ()14j) is negative. We 
want to show that above 7}& the negative eigenvalues of the matrix (|14jl are of higher order 
in 7 and therefore are negligibly small when the coupling is weak. 

Let's consider corrections to the eigenvalues Pi(j+i) and P(i+i)j of the composite state 
f|14jl . The order of magnitude estimate of the smallest one of the corrected eigenvalues is : 
A l i = y/P iU+1) P {i+1)j - jPij/AEij, where AE {j = E^ 1 -E { a + E j b +1 - E{. For simplicity 
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we assume P{{j+i) = P(i+%)j. Then \ l l = 0(P^j + i) — jPij/ AE^-). Below T ib the minimal 
eigenvalue of the state (fTlj) is A_ = 0(— -y/AEu). We shall estimate the ratio r tJ = \ %3 _/\_ 
and show that it is negligible when the coupling is weak. We shall assume without loss of 
generality that the ground state energy is zero: E\ + E\ = 0. Then the partition function 
Z of the composite system is larger than unity It follows that 



x ij 1 p P.,.,,, 7 ZP ZPt^\ 7 p-Eg/T _ -E iU+1) /T 

r lJ = y- = — 3 7 < — = — . (24) 

A - A.E11 AJSu ABu 

We are looking for the maximal value of r %d in the interval < T < T*- 7 , corresponding 
to the condition Piu+i) < jPij/ AEij, i.e. to the negative values of Aj 7 .. T % J is deter- 
mined by the condition A: 7 = 0. The ratio r %3 is positive in the interval < T < T % J* 
and vanishes on its borders. Therefore r lJ has a maximum at < < T*- 7 , which 
is found from the condition dr iJ /dT\ T ij = 0. The calculation gives exp(— AEij/(2T^)) = 
(j/ AEij)(Eij/Ei(j+i)) < (j/AEij) = exp(—AEij/(2T^)), which proves that there is one 
maximum at < < T c u . We remark, that T % c d , corresponding to the largest AE^ over 
all % and j, T* c , is of the order of the upper critial temperature T* c = 0(T UC ). The maximal 
value of r lJ is given by: 



y = AEn f 2E i3 A ^ / 7 \ AE U f 7 A **« 

m 2E ij + AE i] \2E lJ + AE^J \AEtj J ^ 2E ij + AE i} \AE l3 J ' l ; 

where the inequality follows from the fact that 1/e < ( 2 £ +a# ) " < 1 in general. As a 
next step we notice that A£?n < 2i?jj, therefore 

r« < r« < AEn f ^ ^ < A ^ n f ^| ^ (26) 
m 2E ij + AE ij \AE ij ) ~ AE 11 + AE ij \AE ij J 1 ; 

Introducing the definition Xy = AE^ / AEn and taking A£7n = 1, which corresponds to a 
rescaling of the coupling strength 7, leads to: 



_ « . 1 / 7 \ ^ 



Typically the spectrum becomes denser with increasing energy. In that case Xjj = 
AE^ I AEn < 1- Values of corresponding to <C 1 need not be taken into account, 
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because T c lJ < 2]& in this case and as a conseqence X l l > at T > T ih . At = 0(1) the 
upper bound for r scales as 0(7) and therefore corresponding negative eigenvalues of Eq. 
((HJ) are negligible. In this case we expect that « T uc . 

In those cases when x^ = AEij / AEu ^> 1 the upper bound for r scales as 0(1/ Xij) and 
the corresponding negative eigenvalues of Eq. (|14|) can be neglected, too. 

When Xij is moderately larger than unity the upper bound Eq.()27|) for r* J has a local 
maximum. The position of the maximum weakly depends on 7: numerical calculations 
show x^ « 2 — 10 in the range of 10~ 4 < 7 < 1CT 1 . The value of the minimum is a 
monotonically slowly increasing function of 7. In the range 10~ 4 < 7 < 10 _1 numerical 
estimation of Eq. (f2*Tj) shows values 0.04 — 0.1 for the local maximum. It is clear that the 
upper bound Eq. (j27|) for r lJ is far from being tight. In fact, numerical calculations show that 
r y ' is generally much smaller. As a consequence, the corresponding negative eigenvalues of 
Eq. dHJ) can be neglected. 

It can be argued that although each one of the negative eigenvalues of Eq. (|T4"|) is negligible 
at T > the (free) entanglement of the state cannot be neglected. In fact, the minimal 
negative eigenvalue of the partially transposed matrix is not a measure of entanglement. 
Various measures of entanglement have been defined [3(|. In the present context we will 
employ a quantitative counterpart of the Peres-Horodeckii criterion, the negativity 27 1: 



where 



X 



N(p(t)) - mTa J l - , (28) 
TrVX^X is the trace norm of an operator X. The negativity of the state 



equals the absolute value of the sum of the negative eigenvalues of the partially transposed 
state. When the negativity of a composite bipartite state vanishes there is no free entan- 
glement in the state. It can be shown by the order of magnitude analysis similar to the 
analysis above that values of the negativity of the composite state, corresponding to the 
partial transpose (|T4*|) . are generally dominated by the minimal negative eigenvalue. As a 
consequence, the negativity of the state, evolving from the initial thermal product state at 
the temperature T > Tib, is negligible under the weak coupling assumption. 

Figures El and |U display results of numerical calculations of the time averaged negativity 
of the composite state (jOJ as a function of initial temperature for two different kinds of 
unperturbed spectra of the composite system A® B. Both A and B are four level systems. 
The composite system evolves from the initial product of thermal states of A and B under 
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FIG. 3: The time averaged negativity as a function of initial temperature. The composite system 
is constructed from two interacting four level subsystems. The initial state is a product of thermal 
states. The evolution is generated numerically by the Hamiltonian (J2J) (for details of the Hamil- 
tonian see the text) with 7 = 0.05. The dashed lines correspond to the lower bound temperature 
Tib, Eq. (|22|) . and to the upper critical temperature T uc , found numerically. It can be seen that 
the entanglement is vanishingly small in the interval 7}& < T < T uc . 

the Hamiltonian 

Fig. El presents the results of calculations for the following choice of the unperturbed 
spectra of H a and H 6 : e] 1 ' 2 ' 3 ' 41 = {1,5,8,10} and E\ = s/W a . Care was taken to avoid 
resonances and the spectra were chosen to become denser with increasing energy. The 
interaction terms in the Hamiltonian were restricted to (V a ,b)ij — ^»j±i an d the coupling 
strength 7 = 0.05. We see that T uc « and the time averaged negativity (N(p(t))) is 
negligible in the interval Tj& < T < T uc as expected. 

Fig. HJdisplays the time averaged negativity (N(p(t))) as a function of initial temperature 
of the composite state of two interacting four level subsystems A and B with the unperturbed 
energy spectra i?] 1 ' 2 ' 3 ' 4 ^ = {1,3, 7, 13} and E\ = \fE\. The composite state evolves from the 
initial product of two thermal states under the Hamiltonian (J2J), where (V a ,6)y = $i,j±i an d 
the coupling strength 7 = 0.05. In choosing the unperturbed spectra care was taken to avoid 
resonances and to ensure that the maximal value of Xij = AEij/AEn equals the position 
of the local maximum of the upper bound (J27J), corresponding to 7 = 0.05. Fig. El shows 
that the time averaged negativity (N(p(t))) is negligible in the interval 7]& < T < T uc as 
expected. The value of = T* c (the definition of is given after Eq.(J2U)), corresponding 
to the maximal value AE max = max^AE^) is calculated. T* c is in good correspondence 
with the value T uc , calculated numerically. 
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FIG. 4: The time averaged negativity as a function of initial temperature of the composite system. 
The composite system is constructed from two interacting four level subsystems. The initial state 
is a product of thermal states. The evolution is generated numerically by the Hamiltonian (J2J) (for 
details of the Hamiltonian see the text) with 7 = 0.05. The dashed lines correspond to the lower 
bound temperature Tn> Eq. (|22|). to the numerical value of the upper critical temperature T uc and 
to the value T* c , corresponding to the largest spectrum spacing AE max . We see that entanglement 
is vanishingly small at Tu, < T < T uc , as expected, and T* c is a good approximation to the upper 
critical temperature T uc . 

III. ENTANGLEMENT BETWEEN TWO NONINTERACTING SYSTEMS IN 
CONTACT WITH A COMMON THIRD PARTY 

The dynamics studied is of the composite system A® B ®C where systems B and C do 
not interact directly (Cf. FigEJ). The entanglement explored is of the reduced composite 
system B ® C. 

The evolution is generated by the following Hamiltonian: 

H totaZ = H + T V, (29) 

where H = H a + H;, + H c and V = 7V a ® (V& (g> 1 + 1 ® V c ). The analysis is carried out 
in the interaction picture. The initial state is taken to be p(0) = p a ® p b ® p c , where p a , 
p b and p c are thermal states. Since B and C are noninteracting entanglement will appear 
only in the second order in the coupling. Up to second order in 7 the state of the composite 
system A <g> B ® C becomes: 
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Entanglement ? 



FIG. 5: Scheme of interaction for two noninteracting systems in contact with a common third 
party. 



where 



Jo L 



- 7 



o ./o 



v(0, v(O,p(0)' 



dt'cif", 



V(t) = e !ffi Ve 



(30) 



(31) 



In what follows the tag above the p(t) is omitted. 

Next the system is reduced to B®C by taking the partial trace of p(t) over the system A 
degrees of freedom and the partial transposition with respect to the subsystem B is taken: 



p£(0 = p£(0) + M(*), 



(32) 
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where p bc (t) = Tr a (p(t)) and 



M 



t 

->', I Tr a 
'o 

ft rf 

r / / Tr a 

'o Jo 



V(f),p(0) 



dt' 



V(t'), V(t"),p(0) 



(33) 



Let = |z) (g> | A;) be the local orthonormal basis of the system B ®C composed of the 
eigenstates of the Hamiltonian + H c : 



H b , c \i) = E* \i) 



(34) 



where E l bc , i = 1,2, is the unperturbed energy spectrum of the Hamiltonian Hf, c . Since 
PbM = Pb® Pc- 



PbM Tb \ik) = p b M \ik) = P b ®Pc \ik) = P ik \ik) 



(35) 



where P ik = Pb,iPc,k, and p a ,i,Pb,k are defined by p b)i = (i\p a \i) and p C)k = (k\p b \k). The 
matrix elements of p bc (t) Tb are given by: 



(ik \p bc (t) Tb \ jl) = Pi k S(ik),(ji) + M t 



ik,jl ; 



(36) 



where by definition M, 



ik,jl 



ik 



M 



From this point the calculations proceed along the same lines as in Section II following 
Eq.ffTT]). The minimal eigenvalue of the partially transposed reduced state p bc (t) Tb is shown 
to be negative at sufficiently low temperatures and the lower bound temperature T lb is 
calculated. 

The negative eigenvalue of the partially transposed composite state Eq. (J32j) is calculated 



to the leading order in the coupling strength 7 assuming in 



V, 

Section II the eigenvalue is found from the spectrum of the 2x2 matrix: 

P\2 + Mi2,12 M\2,2l 

-^12,21 °21 + ^21,21 

completely analogous to the matrix (|16|). The eigenvalues of Eq. ()37|) are: 
P12 + M 12 12 + P21 + M2121 



mi) oc 5 ni m t ±i- As in 



(37) 



A- 



(38) 



± 



y/(P 12 + M 12 ,12 + P21 + M 21 ,2i) 2 - 4((P 12 + M 12 , 12 )(P 21 + M 21 , 21 ) - |M 12 ,2l| 
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and the eigenvalue A_ becomes negative when (Pi 2 + Mi2,i2)(-f2i + ^21,21) < |^i2,2i| 2 - 

To calculate Mj.2,12, ^21,21 and Mi 2 ,2i we first note that the integrand in the first order 
term in Eq. is: 



v(0,p(o) 



V b (f),Ph 



V 6c (t'),P 6c (0)_ 

V c (t'),Pc 



Pb 



(39) 



V, 



v fe (t') T ,p b 



Pfc 



v c (t'),p c 



where \ V a j means the thermal average of the operator V a and the notation Vf, c = 1 + 
1 ® V c is introduced. The initial condition p bc (0) = p b (g) p c was used. Since p b c |i) = 5^1 |i) 
the term Eq. (J39)) does not contribute to the eigenvalues of the matrix (J37)) in the first order. 

To simplify the calculation of the second order corrections it is assumed that the thermal 
average of the system A coupling operator ( V ) vanishes. This assumption is not crucial 



for the qualitative picture of temperature dependence of the entanglement. Moreover, it 
is in line with common models of coupling, for example, the Caldeira-Leggett model |31[, 
dipole interaction with the electromagnetic field (3^], etc. The integrand in the second order 
term in Eq. (jHB^l is: 



Tr a 



V(f), V(t"),p(0) 



v o (0v a (O) V bc (t') } V bc (t")p bc 



v a (t")v a (0> v bc (t'),p bc v bc (t") 



(40) 



Expanding the thermal averages in the orthonormal basis |n) of the Hamiltonian H a leads 
to: 



V(f), V(t"),p(0) 



m,n 



rn 



V,. 



n 



(cosK^' - O) v fc (f), v te (0, 



+ 2sin« n (t'-t")) V 6c (t'),{V 6c (*"),p 6 , 



where is the energy difference between the states \n) and |m) of the Hamiltonian H a , 
the jx, Y j designates anticommutator of operators X and Y and p a ^ n = (p a )„„. 

For simplicity the notation C(t', t") is used for the operator (JJTJ). Expressing the operator 
Vfc c in terms of Y b and V c we put the matrix elements of C(t',t") into the following form: 
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(12|C(t',t")|12> = -2Pn^pa, n 

m,ra 

(21|C(t',t")|21> = -2P n ^ Pa , n 

m,n 

(12|C(t',t")|21> = PiiJ2p^\ 

m,n 
-i(uj b t"+ui c t')^ 



III 



m 



V, 



n 



n 



n 



V, 



co${{u a mn + u c ){t'-t")) 

cos(« n + a; b )(t'-i")) (42) 
\ A v c 1^ e ' wL( ''" t " ) (e" i(uti,+Ut '" ) 



+ e" 



where a;& )C stands for the energy difference between the first excited and the ground states 
of the unperturbed subsystem B ( C). The matrix elements Mi 2 ,i2, M 2 i t 2i and Mi 2j2 i are 
given by: 



M 12 , 12 = - 7 2 / / (12| C(f,t") \12) dt'dt". 
Jo Jo 

M 2 i, 2 i = - 7 2 / / (21|C(t',t")|21) dt'dt" 
-t /.tf 



(43) 



M 



12,21 



"7' 



/ / (12\ C(t',t")\21) dt'dt". 
Jo Jo 



The integration is straightforward but the final expressions are cumbersome. Two cases 
are considered explicitly: (a) >> Ub tC and (b) Ub tC » ^mn- m both cases it is shown 
that at sufficiently low initial temperature of the system B®C one of the eigenvalues of the 
matrix (|3*7)l is negative and the lower bound temperature T Lb is calculated. 

A. Two "slow" systems interacting with a "fast" common third party 

Performing the integrations in Eq. (}43|) and taking the leading terms in ujb )C /^mn brings 

to: 



M 



12,12 



M> 



21,21 



M 



12,21 



m,n 

m,n 
2l 2 ^2p a ,n 



m 



m 



m 



V, 



V, 



n 



n 



n 







2 > 







2 > 



2sin(« n + o; c )t/2) 1 



a 



2 sin(« n + ^)t/2)^ 



a ^2 



(44) 



12 



f\ _ e —i{^b+w c )t\ 
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At T 



the minimal eigenvalue of Eq. 



is given by A_ = 
--^■^12,12-^21,21 - |^i2,2i I 2 , which to the leading order in uj^Juj^ gives A_ = -| ^12,21 1 2 ■ 
This proves that the system B g) C becomes entangled at sufficiently low temperature. We 
note that this result holds at any finite temperature of the system A. At infinite temperature 
of the system A M\ip\ = and no free entanglement is generated in the system B ®C. 

At finite initial temperature of B <g> C the condition A_ < translates to P\iP%\ < 
7 4 |M 12] 2i| 2 -Pii to the leading order in uj^^l^mn- The lower bound temperature T ib is found 
from the condition P12-P21 = 7 4 1 ^12,21 \ 2 P\i- Since |Mi2,2i| is an oscillating function of time 
the amplitude of |M 12j 2i| must be substituted for |M 12) 2i| in this equality, which leads to the 
following equation defining the lower bound temperature: 

; 2 



4 7 2 



V, 



exp( 



u b + UJ C 

U! b + UJ c 

2T, h h 



E 



Pa,n 



III 



n 



'P2P3 



, ,0 



finally leading to: 



T, 



-{u h + u c ) 



lb 



2 In 4 7 



(2 V b l)(2 V e l) 



E 



(rn V a nj 



m,n 



(45) 



(46) 



A generalization of the formula to the case of interaction of the form ^ 7; V a ® (V 6 <8> 1 + 
1 ® V c ) can be carried out along the same lines. 

The entanglement in the reduced system of two noninteracting "slow" spins interacting 
with the "fast" four level "bath" was explored numerically and the results are plotted on 
Fig. El The shaded area in the parametric space of the logarithm of inverse coupling 
strength and the inverse initial temperature of the spins represents parametric values for 
which no entanglement develops in the course of the evolution. The border of the shaded area 
corresponds to the critical temperature for various coupling magnitudes. The Hamiltonian 
of the composite system is: 



H = H a ® l b ® i c + ^(i a (8) (& b z <g> i c + V21 b <8> &D) 

+ 7 v a ® {& b x ®i c + i b ®& c x ), 



(47) 



where (H a )^ = E^' 2 ' 3A} = {0, lOu, 20u, 30u} and (V a ) y = The temperature of 

the thermal initial state of the "bath" is T = 5u. The value of u chosen for the numerical 
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FIG. 6: The shaded area in the parameter space of the inverse initial temperature T of the "slow" 
spins and the logarithm of the inverse coupling strength 7, represents values of T and 7 where 
entanglement does not develop in the course of the evolution. The composite system of two "slow" 
spins interacting with a " fast" four level system evolves from the initial product of thermal states 
under the Hamiltonian (|47|). The dashed line is the plot of 2)&, Eq. (jl6|) . Up to the coupling 7 = 1 
its correspondence to the border of the shaded area is very good. 

calculation is unity. The correspondence of Eq. (J46|) (the dashed line) to the numerical 
values is very good up to a coupling strength of the order of unity. We note that for large 
values of the coupling strength 7 the critical temperature asymptotically tends to a finite 
constant value. 

B. Two "fast" systems interacting with a "slow" common third party 

The case u^ c » o;^ n is more complex. To demonstrate entanglement at zero temper- 
ature of the system B <g> C two simplifying assumptions were added. The first is that the 
temperature of the system A is also zero. The second is that the matrix elements of V a 
couple only the neighboring states: (n| V a \m) oc <5 n ,m±i- Under these two assumptions the 
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expressions for Mi.2,12, ^21,21 and Mi 2< i 2 Mm 21 2i — |^i2,2i| 2 become: 



M 



12,12 — P11 



2 7 



21,21 



P11 



2 7 



sin 



(jjjg+u c )t 



V, 



sin 



M 12jl2 M M2121 - |M 12j21 | 2 



p2 

-'11 



^ 2 7 2 



LO b UJ c 



V, 



(48) 



where 



S'(t) = sin(u;at)[sin(u;bt) + sin(u; c t) 
- sba.((u a + u b + u) c )t)]. 



(49) 



To estimate S(t) new variables x = sin(u a t), y = sin (u;&t) and z = sin(u; c t) are introduced. 
Ignoring the zero measure set of commensurable frequencies we can treat the function S(t) 
as function of three independent variables x, y and z. The range of S(t) in the cube, defined 
by — 1 < x, y,z < 1, can be explored numerically and is found to be: s < S(t) < 3, where 
s « -1.6834. Therefore, from En. pBj) M 12jl2 M M2121 - |M 12j21 | 2 < 0, which proves that at 
zero temperature A_ < (Cf. Eq.(|55jl) and the systems B and C are entangled by the 
interaction with the system A. 

The lower bound temperature is determined by the condition A_ = 0, which translates 
to (P i2 + Mi 2j i 2 )(P 21 + M 2 i j21 ) = \M 12 ,2i\ 2 (Cf. Eq. pSjl ). The latter condition can be 
put in the form (M 12il2 M 21 , 21 - \M 12i21 \ 2 ) + P 12 P 21 + Pi 2 M 21 , 21 + P 21 M 12 , 12 = . Since 
M 12 21 and M 2121 are nonnegative independent functions of time the minimum value of 
(M 12i12 M 21i21 -|M 12>21 | 2 )+P 12 P 21 + P 12 M 21i21 + P 21 M 12i12 is obtained at M 21>21 = M 12jl2 = 0. 
Then the lower bound temperature can be calculated from the condition that the amplitude 
of Mi 2) i 2 M 2 i, 2 i - |ilPi2,2i| 2 equals -Pi 2 P 2 i: 

(l V c 2^ (l V b 2) f^P~ 

(50) 



LU b UJ c 



12-r 21 
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finally leading to: 

l 'VII OJ b Wc I 

It is interesting to note that in this case does not depend on the time scales of the "slow" 
system. 

The entanglement in the reduced system of two noninteracting "fast" spins interacting 
with the "slow" four level "bath" was explored numerically and the results are plotted on 
Fig. [3 The shaded area in the parametric space of the logarithm of inverse coupling 
strength and the inverse initial temperature of the spins represents parametric values for 
which no entanglement develops in the course of the evolution. The border of the shaded area 
corresponds to the critical temperature for various coupling magnitudes. The Hamiltonian is 
chosen to be similar to the previous example, Cf. Eq. (j47|) . but time scales of the subsystems 
are reversed: 

H = H a <g> l b <g> i c + 5a;(i a <g> (& b z ®l c + V21 b <g> &D) (52) 

+ 7 v a ® (& b x ®i c + i b ®& x ), 

where (H. a )ij — dijEai Ei 1 ' 2 ' 3 '^ = {0, oj, 2u, 3a;} and (V )jj = 5^. The temperature of the 
thermal initial state of the "bath" was chosen as T = 0.01a;, which is small compared to the 
energy scale of the "bath" chosen for the numerical calculation: to = 1. The dashed line on 
the Fig. is a plot of Eq. (joTj) and the correspondence to the border of the shaded area at 
coupling strength up to the order of unity is good. 



IV. SUMMARY AN CONCLUSIONS 



Entanglement is created by both direct and indirect weak interaction between two initially 
disentangled systems prepared in thermal states at sufficiently low temperatures. The study 
is restricted to the conditions where the ground states of both systems are not invariant 
under the interaction and the interaction is nonresonant. As a consequence, the present 
analysis left out some interesting models such as the Jaynes-Cummings model 33j. The 
Jaynes-Cummings model of interacting two level system and a quantized field mode was 
investigated in Ref.[16]. It was found that no free entanglement is generated in the course 
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FIG. 7: The shaded area in the parameter space of the inverse initial temperature T of the "fast" 
spins and the logarithm of the inverse coupling strength 7, represents values of T and 7 where 
entanglement does not develop in the course of the evolution. The composite system of two "fast" 
spins interacting with the "slow" four level system evolves from the initial product of thermal states 
under the Hamiltonian (|52|) . The dashed line is the plot of 2)&, Eq. (|51j). Up to the coupling 7 = 1 
its correspondence to the border of the shaded area is good. 



of the evolution of the composite system if the initial temperature of both the subsystems 
is sufficiently high. 

The generation of entanglement in cases of the weak resonant direct and undirect inter- 
actions will be treated separately 3^ |. 

In the case of indirect interaction to show entanglement at T = we have assumed that 
the thermal average of the third party coupling term in the initial state vanishes. The 
reason for the assumption was technical. It should be noted that many system-bath models 
of linear coupling satisfy this assumption j^j] • The additional technical assumption was that 
the coupling terms of the noninteracting parties possess matrix elements only between the 
adjacent energy states. Here, too, the assumption is general for weak coupling models. Two 
cases of time scale separation were considered explicitly. The first is the case of two "slow" 
systems interacting via the "fast" third common party. The second is the case of two "fast" 
systems interacting via the "slow" third common party. In the first case the entanglement 
was shown to appear at sufficiently low initial temperature of the "slow" systems for any 
finite temperature of the third party. In the second case the entanglement develops at 
sufficiently low initial temperature of the "fast" systems. In this case we assumed that the 
third party was prepared at zero temperature and that the third party coupling agent has 
nonvanishing matrix elements only between the adjacent energy states. This assumption is 
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stronger than just assuming that its thermal average vanishes. 

In these cases of indirect interaction and in the case of the direct interaction between the 
parts we have shown that if the initial temperature of the bipartite state is zero entanglement 
is generated by the interaction. At sufficiently high temperature the composite state remains 
PPT in the course of evolution. From these results it follows that a lower critical temperature 
Tic exists: if the initial temperature of both thermal states is below T\ c the interaction 
generates entanglement in the course of the evolution, and if the initial temperature is 
sufficiently close to T; c from above the the composite state remains PPT forever. When the 
composite system is finite dimensional there exists an upper critical temperature T uc : if the 
initial temperature of both thermal states is higher than T uc the composite state remains 
PPT in the course of evolution and if the initial temperature is sufficiently close to T uc from 
below entanglement is generated. We conjecture on the basis of numerical experiments that 
Tic — T uc m general. In both cases of a direct and an indirect interaction between the initially 
disentangled systems, prepared in thermal states, we calculated the lower bound Tj& for the 
lower critical temperature TJ C . When the initial temperature of both thermal states is below 
Tb the interaction generates entanglement in the course of the evolution. For temperatures 
above the lower bound T}& the negativity of the partially transposed composite state is zero 
in the leading order in the coupling strength and therefore negligible in the weak coupling 
limit. It follows, that may be considered as the physical critical temperature for the 
negativity of the composite state. 

Separable states can be considered as classical states, because they lack quantum cor- 
relations. One may hope that, as a consequence, the dynamics of separable states can be 
efficiently simulated on classical computers. Whether this is possible is an open question 
in quantum information science. If a moderate scaling procedure exists for the simulation 
of the dynamics of a separable bipartite state, then it seems that such a procedure exists 
also if the evolving state remains PPT for all times. Ref.j^ has proved that a density 
operator p supported on a M x N dimensional Hilbert space (M < N) with positive partial 
transpose and a rank smaller than or equal to N is separable. It follows that a PPT state 
of dimension N is always separable when embedded in the larger Hilbert space of dimension 
iV 2 or higher. The dynamics of the low dimensional PPT state will be physically equivalent 
to the dynamics of the high dimensional separable state which can (hopefully) be simulated 
efficiently on the classical computer. 
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The present analysis shows that above a critical temperature T} 6 the PPT character of a 
composite state is preserved along the evolution. The challenge is to construct an effective 
simulation for the dynamics of a composite quantum systems at finite temperature employing 
classically based computers. 
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